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We propose a method for acting on the spin state of a spin-| localized particle, or qubit, by 
means of a magnetic signal effectively generated by the nearby transit of a magnetic soliton, there 
conveyed through a transmission line. We first introduce the specific magnetic soliton of which we 
will make use, and briefly review the properties that make it apt to represent a signal. We then 
show that a Heisenberg spin chain can serve as transmission line, and propose a method for injecting 
a soliton into the chain by acting just on one of its ends. We finally demonstrate that the resulting 
magnetic pulse can indeed cause, just passing by the spin-| localized particle embodying the qubit, 
a permanent change in its spin state, thus realizing the possibility of getting through to a single, 
localized qubit, and manipulate its state. A thorough analysis of how the overall dynamical system 
operates depending on the setting of its parameters demonstrates that fine tuning is not necessary 
as it exists an extended region in the parameters space that corresponds to an effective functioning. 
Moreover, we show that possible noise on the transmission line does not invalidate the scheme. 


I. INTRODUCTION 

The ability of addressing, initializing, and possibly con¬ 
trolling one single qubit without spoiling its quantum 
features or disturbing other nearby qubits is a necessary 
prerequisite for putting a quantum device into operation. 
Depending on the specific device architecture, however, 
this can be a most challenging task, as it implies the 
opening of the qubit towards an environment that, in a 
way or another, embodies some macroscopic apparatus. 
One possibility for avoiding that this opening alter the 
fragile properties of the qubit is that of placing the ap¬ 
paratus at a distance, and use a transmission line for 
conveying a proper signal to the qubit itself. In particu¬ 
lar, when the qubit is represented by a localized magnetic 
particl^i®^ it comes quite natural that the above signal 
be a time-dependent magnetic field, which nonetheless 
leaves the question open on how to realize the transmis¬ 
sion line. To this purpose, we here propose the use of 
classical spin chains featuring soliton-like excitations, a 
choice suggested by these observations: i) A soliton faith¬ 
fully represents a signal in so far as it is a finite-energy 
excitation which is well localized in space at any given 
time, and can travel at fixed velocity with constant pro¬ 
file; ii) Solitons are known to travel undisturbed along 
their medium, which allows us to put the apparatus that 
generates the pulse at great distance from the qubit and 
yet be sure that the signal will pass near its target un¬ 
deformed; in) Solitons relative to the very same model 
can have different shapes and energies, which gives us 
the freedom of choosing the signal that best controls the 
qubit, without modifying the transmission line. 

As for this latter component, we know that some classi¬ 
cal helds dehned on a one-dimensional space display soli- 
tonic excitations, whose renowned stability stems from 
the competition between linear and non-linear terms in 
the held’s equations of motion (EoM). Based on the well 
established connection between classical vector-held the¬ 
ories and models of interacti^ spin-S' particles on dis¬ 
crete one-dimensional latticed, one can expect soliton- 
like excitations to typify some spin chains^fUl, which di¬ 


rectly suggests that one such chain can serve for support¬ 
ing a solitonic signal. Moreover, as extensively shown in 
the literatur^i^HIZl ^ renormalized classical approadpsHUl 
is often appropriate for describiM the actual behaviour 
of real compounds with S' > 1 which allows us to 

treat the signal’s generation and propagation at a classi¬ 
cal level. 

The paper is structured as follows: In Sec. |lT] we in¬ 
troduce and characterize the soliton of which we will 
make use as non-linear excitation of a classical vector 
held. Based on the fact that such held is the contin¬ 
uum limit of a spin model on a one-dimensional lattice, 
namely the Heisenberg spin chain, in Sec. |III| we numeri¬ 
cally check that dynamical conhgurations corresponding 
to discrete versions of the above solitons can be solu¬ 
tions of the chain’s EoM; in the same section we pro¬ 
pose a method for injecting solitons by applying a time- 
dependent magnetic held to one end of the chain. In 
Sec. |IV| we thoroughly analyze the effects of the soliton’s 
transit on the spin state of the spin-^ localized particle 
embodying the qubit, considering both an ideal and an 
injected soliton, as well as the possibility of thermal noise 
along the chain. Comments upon possible experimental 
implementations of the scheme are put forward in Sec. 
with attention focused on the validity of the assumptions 
we have made in order to make our scheme function. 


II. SOLITONS AS SIGNALS 

Consider the classical vector-held in one spatial di¬ 
mension, s{x, t), such that |s(a:, t)| = 1, with Hamiltonian 
density 

n{x) = \jS^[d,,s{x,t)Y+-tHS[l-s^{x,t)] : (1) 

its EoM 

dts{x, t) = s{x, t) X [JS d^s{x, t) -I- yff] , (2) 

with H = ( 0 , 0 , 77 ), have been showrPS! to admit ana¬ 
lytical solutions, corresponding to localized excitations. 


2 



FIG. 1. for tan/3 = 2: = 1—2sin^/3/ cosh^^ as a 

function of x—vt and (inset) the corresponding trajectory of 
the in-plane components and . 

stable under collision^^, that travel at constant velocity. 
In polar coordinates these solutions read 


dependence of the spin-variables will be hereafter under¬ 
stood, whenever possible. The Heisenberg chain is de¬ 
fined by the Hamiltonian 

Hchain = -jJ2 '^'■'^'+1 - 

I I 

JQ2 

= + jHS ^(1-sf) + const , 

I I 

(9) 

where J > 0 is a ferromagnetic nearest-neighbor cou¬ 
pling, H = (0, 0, H) is an external magnetic field, and 
7 is the gyromagnetic ratio. In the (continuum) limit of 
vanishing lattice spacing d—^0 with finite S/d and Jd^, 
the model § reproduces Eq. 0 with 

S = ^ and J = Jd^ . (10) 

d 

The EoM for si, obtained from the Poisson brackets 
{sf,Sj} — S~^ Sij 6°"^^ s/, consistently have the same 
form of Eqs. ([^, reading 

dtSi = SiX [JS{si+i+Si_i) + -fH] . (11) 


where 


6^^\x,t) = 2 sin ^ (sin/3 sech^) , 
ip^^\x,t) = v^o + cot/3 ^ -I- tan“^ (tan/3 tanh^) 

(3) 

^ x—vt X t ... 

e = ^ , (4) 

A/3 Tp 


and the parameter /3 univocally characterizes each u > 0 
soliton, setting its characteristic 


amplitude : 

V 

2B = 2 arccos —, , 

2^JSyH 

(5) 

length : 

- {^H sin/3 ’ 

(6) 

energy : 

ep = %S\/ JSyH sin/3 , 

(7) 

time : 

1 

~ 7Lfsin2/3 ' 

(8) 


A dynamical {v > 0) soliton defined by Eqs. ^ will be 
hereafter referred to as ”/3-soliton”. Notice that Eq. ([^ 
sets a maximum value for the velocity, |u| < 2^ J S')!!, 
implying that the second term of the Hamiltonian Q 
must be finite in order for the model to support dynam¬ 
ical solitons. Once this condition is fulfilled, a /3-soliton 
can be readily seen as a signal, i.e., a field’s dynamical 
configuration with a distinctive trait that can be spotted, 
for time intervals of the order of Tp, in a spatial region 
of size \p that moves with constant velocity in the one¬ 
dimensional space where the field is defined. An example 
of /3-soliton is shown in Eig. 


III. HEISENBERG CHAIN AS TRANSMISSION 
LINE 

Let us now consider a classical spin chain, i.e., a one¬ 
dimensional array of (spin-)vectors Si = S si, whose 
magnitude S has the dimension of an action; the time 


Despite the analogy, analytical soliton-like solutions of 
Eqs. (Ill are not known; however, as the continuum ap¬ 


proximation does make sense whenever the relevant con¬ 
figurations vary slowly on the scale of the lattice spacing, 
we expect that, for Xp ^ d, the discrete counterpart of 
a /3-soliton, s[^^(t) defined by Eqs. 0 with 


C = l\l sin/3 — 7iLtsin2/3 , 

J o 


( 12 ) 


might still represent an excitation of the model ([^. In 
fact, by numerically solving Eqs. (Ill, we have checked 


that the Heisenberg spin chain properly supports /3- 
solitons whenever the Zeeman energy ^HS is much 
smaller than the bond energy JS"^, as implied hy Xp':S> d 
via Eqs. 0 and ( [To| . This result fits with the exper¬ 
imental observation, in quasi one-dimensional systems, 
of magnetic behaviours whose origin can be unequivo¬ 
cally ascribed to the presence of soliton-like excitation^^. 
In what follows, we will therefore take that the time- 
dependent chain configuration {s/(f) = s|^^(t)} is a so¬ 
lution of the discrete EoM embodying the signal we 
want to convey to the qubit, with the respective Heisen¬ 
berg chain serving as transmission line. 


A. Injecting a soliton into the chain 

Let us now consider the problem of making one specific 
soliton to exist and run through the Heisenberg 

chain. In the process of accomplishing this goal, we first 
notice the following: Consider a finite (though long at 
will) chain, with (2L -|- 1) spins sitting on sites labelled 
from —L to +L. Suppose S-L{t) evolved as if a soli¬ 
ton were reaching it travelling from a fictitious, infinitely 
left-extended chain, I < —L: that soliton would continue 
travelling towards the region I > —L, i.e., it would be 
successfully injected into the chain, at least in the con¬ 
tinuum limit. Therefore, enforcing 

s.ut) = sL^l(t) 


(13) 
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as a boundary condition should result in the selection of 
the configuration corresponding to (t) , amongst all 
those that solve Eqs. (111. On the other hand it can be 
easily seen that Eqs. (11) with condition (13) enforced 
are the EoM of a Heisenberg chain with —L < I < L, 
and an auxiliary magnetic field 


IV. QUBIT DYNAMICS 

A qubit is a physical system that can be described by 
the spin-i operator represented by the Pauli matrices 
& = (7^, (T^); its state p{t) in terms of the Bloch 

vector n{t) = Tr[/3(t) < t ] , reads 


(14) 


Kt) = \ + ■ 


(15) 


acting just on S-l+i- This suggests that, by applying 
the magnetic field (14) to one end of the transmission 
line, it should be possible to generate a soliton-like sig¬ 
nal, that will then travel through to its target. The con¬ 
sistency of the above description has been checked, for 
different values of /3, as follows. 

First, we have numerically solved Eqs. 0 with the 
time-dependent constraint Eq. ( 0 ), and 2Ld^Xp, by 
means of a second-order symplectic algorit hnpIHSU The 
chain has been initialized in the ferromagnetic config¬ 
uration, {si=z}, as well as in some possible ther¬ 
mal configurations. These have been determined, re¬ 
ferring to the quadratic approximation of the Hamilto¬ 
nian 0 which is diagonal in Fourier space with fre¬ 
quencies ujk = 2JS'(1—cosfc) -|- yH, as inverse Fourier 
transforms of generated sets of independent variables 
{Sfe.Sfe} with variances (s^s^fe) = (SfeS^.^) = k^T/{SuJk). 
The resulting configurations have the thermal correla¬ 
tors ((sf-sf+i) 2 ) = ~ r and ((sf)^) = 

((sf)2) ~ T^JShH, where T = k^T/JS'^. 

We have then analyzed the resulting solutions and 
found that when the field is applied to 


i.e., after the injection of dynamical configura¬ 

tions r;(/3;t), identifiable as soliton-like, actually appear 
into the chain. In order to better characterize these con¬ 
figurations, we have numerically measured their veloc¬ 
ity v' and, assuming the validity of Eq. ([^, we have 
obtained values for the respective amplitudes, 2j3' = 
2 arccos[f'/(J5'7H)]. These values have been found to 
agree with those independently determined by fitting 
with Eqs. for all values of /3 considered. 
Moreover, by monitoring the chain’s energy throughout 
the numerical integration, we have calculated the total 
work made by the forcing term, and found it very close 
to ep', meaning that the work done on the chain does 
actually correspond to the energy needed to generate a 
soliton ^ (t) . Summarizing, the above analysis con¬ 

firms that: 


by applying the field 6 *^^^ (t) defined in Eq. (14) to 
the left tail of the chain, 


• a dynamical configuration is generated in¬ 

side the chain itself, 

• with the essential features of a soliton Sj ' {t). 

The above picture is also confirmed for T > 0, and even 
for rather narrow injected solitons (A^ ~ 5d). In Fig. 
we show the colour-density plot of Tf (/3; t) as a function 
of I and JSt, for different T. The strong resilience of the 
generated signal is evident: even when fully embedded 
and barely recognizable within the thermal noise, as in 
the last panel, its time propagation along the chain can 
still be easily followed. 


In our scheme, the qubit is realized by a localized spin- 4 
particle, sitting near the site of the chain labelled by the 
index “0”. We assume that the way the qubit feels the 
presence of the chain’s magnetic moments, as depicted in 
Fig.i can be generally described as a Zeeman interaction 
with an effective magnetic field proportional to 


s(i) = , ( 16 ) 

3 

where Sj are the unit vectors entering Eq. 0, and pj 
is expected to decrease rapidly with |j|. In fact, the de¬ 
tailed dependence of pj on j is not relevant, particularly 
if, as in the present scheme, the time dependence of the 
magnetic moments is primarily due to the transit of a sig¬ 
nal whose length is of a finite number of lattice spacings. 
Therefore, we can safely choose a Gaussian dependence 
Pj = Aexp{—j^/2a'^) , where A is such that J2jPj = 1 j 
and the standard deviation a characterizes the interac¬ 
tion range, in units of d. 

From Secs. and |HI| we have learned that the pres¬ 
ence of a constant and homogeneous magnetic field H is 
necessary for the Heisenberg spin chain to support soli¬ 
tons with finite velocity: therefore, we take H 0 and 
identify its direction with the quantization axis used for 
encoding the qubit states into the spin degree of freedom 
of the spin-1 localized particle. The qubit’s Hamiltonian 
thus reads 

Tiquhit = -[laH z +gS s{t)] , (17) 

where ja is the gyromagnetic ratio of the particle real¬ 
izing the qubit, and g is an overall coupling constant. 
The corresponding evolution of the qubit’s Bloch vector 
is ruled by the equation 


OtTI = n X [S z + p s(r)] , (18) 

where z = (0,0,1), and r = jHt is the (chain’s) di¬ 
mensionless time that will be hereafter used. The two 
dimensionless parameters 




(19) 


characterize the qubit’s interactions with the external 
and the effective field, H and s{t), respectively. No¬ 
tice that, despite the chain parameter 7 does not appear 
in the qubit’s Hamiltonian ( |T7| , it does enter the EoM 
for the qubit’s Bloch vector via the definition of the di¬ 
mensionless time r; in fact, the relevant time scale of the 
overall dynamics is exclusively set by the chain Hamilto¬ 
nian 0 , a statement based on the implicit assumption 
that the presence of the qubit has no effect (no ‘back- 
action’) on the chain itself. We will further comment 
upon this assumption in Sec|V| 
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FIG. 2. Samples of generated soliton-like excitations in a discrete chain of 500 spins; the parameters of the injected 

soliton are 'yH/JS = 0.05 and tan/3 = 2. Density plots are shown for the space-time evolution of sf{t) at zero (upper 

left panel) and finite temperature, indicated in the lower, right edge of each panel. The propagating soliton is reflected by the 
open boundary at site n = 500. The thin, red line reports the time dependence of sf{t) at the site n = 250. Thermal noise 
makes the generated soliton-like excitation broader and faster. 
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FIG. 3. The qubit interacts with a bunch of moments of 
the classical spin-chain, with couplings ji = gpi; a constant 
uniform field is applied to the overall system. 


Suppose now that a magnetic signal in the form of 
r/(/3;r) run through the chain (T = 0). In the early 
stage of the process, at a time when the soliton is still far 
from the site 0, it is s{t) (xz and the qubit Bloch vector 
undergoes a uniform precession around z, unless it is not 
initially aligned along the z-axis itself. In order to isolate 
the qubit evolution exclusively due to the soliton transit, 
it is therefore convenient to choose n{Ti)=z, with Ti the 
earlier time when s_]^ yf z. Notice that this does not im¬ 
ply adding a previous single-qubit manipulation step in 
the overall scheme, but rather preparing the whole sys¬ 
tem in a globally aligned state, which is readily obtained 


asH ^0. 

Consider now a time Tf during the final stage of the pro¬ 
cess, i.e. after the soliton has travelled along the chain far 
beyond the qubit: the qubit’s Bloch vector, set in motion 
by the soliton’s transit, can 1) align back to z, 2) tilt- 
up and hence precess around i, or 3) perfectly flip and 
anti-align along —z. Situations 2) and 3) are those in 
which we are most interested, as they represent the pos¬ 
sibility of permanently modifying the qubit state, which 
is in fact the final goal of our scheme. In order to ana¬ 
lyze the conditions under which they are obtained, one 
must numerically integrate Eqs. (181 with the effective 
field s{t) as from Eq. with 


= rj(/3;r) . 


( 20 ) 


For the sake of clarity, in what follows we will specifically 
concentrate on the case when the qubit’s response to the 
signal consists in a permanent flipping. 


A. qubit flipped by an ideal soliton 

Let us first take = sP^(t) and the chain ini¬ 

tially prepared in the ferromagnetic state (T =0). When 
a = 0 we knovi^ that whenever (5 = 0 the qubit always 
goes back to its initial state: therefore, in order to obtain 
a permanent flipping, the physical object embodying the 
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FIG. 4. Contour plots of n^(T{) as a function of S and fi, for a^O. In each panel the magnetic signal acting on the qubit is 
that produced by an ideal /3-soliton, 


qubit must have a finite gyromagnetic ratio. As study¬ 
ing n^(Tf) suffices to distinguish the above situations 1), 
2), and 3), in Fig. [^we plot n^{T{) in the plane 
when S = fi = l the flipping is complete, while the change 
in n^(Tf) decreases monotonically when getting far from 
this point. Remarkably, for 6 = fj,= l there is no depen¬ 
dence on (3: the qubit is flipped whatever the amplitude 
of the signal running through the chain. An additional 
feature, numerically observed and clearly seen in Fig. 
is that n*(Tf) is symmetric in the exchange 6 -H- fj,, even 
though the evolution of the qubit may be different in 
the two cases. The most relevant feature displayed by 
Fig.i however, is the presence of a region where almost 
complete flipping occurs: this means that fine-tuning is 
not necessary and if 6 is difficult to alter one can still 
act on or viceversa, depending on the specific physical 
realization of the scheme. 

When a finite interaction range (a 7 ^ 0) is considered, 
as in the case shown in Fig. the qubit’s dynamics is 
qualitatively similar to that observed for a = 0 ^, but the 
value of n^{Ti) is found quite sensitive to a itself. How¬ 
ever, an almost complete flipping, even better than that 
observed in Fig. can be obtained by further adjusting 
the available parameters. To this respect, notice that the 
ratio h = / JS is a relevant quantity, when a 7 ^ 0 , as it 
contributes to set the length scale of the soliton: For ex¬ 
ample, h = 0.05 and tan/3 = 2,1, 0.5, 0.2 define /3-solitons 


with A /3 = 5d, 6.3d, lOd, 22. 8 d, respectively. 

In Fig. we show a contour-plot relative to n^(Tf) in 
the plane ^an/3,/i), for d = 1 and different values of h. As 
expected, for smaller tan/3 the qubit’s dynamics is less 
affected by the finite interaction range, as broad solitons 
(A /3 » ad) are little modified by the ‘smearing’ entailed 
by Eq. (16). In particular, the plot for h = 0.02 shows 
that the partial flip shown in Fig. [^can be improved by 
taking smaller /S, i.e. longer solitons, or increasing /i, 
i.e. the qubit-chain coupling. The flip quality decreases 
also when, due to the phase term cot/3 appearing in 
Eq. @, the X and y components of s{t) shrink under the 
smearin g (|l^ . This effect can be reduced, as suggested 
by Eq. (|12|), by requiring aVh cosfd, i.e., for small /3, 
h ^ a~‘‘. In fact, Fig.j^shows that by taking a smaller h 
the flip quality can be made to approach optimal values 
in an extended region of the jS-y, plane. Referring to 
the definitions (19), this optimization can be typically 
performed by driving the external field only. 


B. qubit flipped by a generated soliton 

Let us now consider the case when the soliton run¬ 
ning through the chain is not ideal, but rather a gen¬ 
erated one, r/(/3;r) ~ ^(r). In Fig. [^we show the 
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FIG. 5. Time lapse of qubit’s evolution (represented by n(r)) 
while a /3-soliton propagates along the chain. The graphs on 
the left side display s^(r) (full lines) and s^(t) (dashed lines), 
i.e., the components of the effective field acting on the qubit 
as a consequence of the soliton’s transit, drawn up to the same 
time r. The parameter values are: h — 0.01, tan/3 = 0.2, p = 1, 
5 = 1, and a = 3. 


qubit’s state evolution when tan/3 = 2 and = 5 (a = 0): 
in the left panels one can appreciate that the evolution of 
n^(T) follows that of the generated soliton, both for zero 
(top) and finite (bottom) temperature; the right panels 
display the overall trajectory of the qubit’s magnetiza¬ 
tion on the Bloch sphere. The qubit’s behaviour under 
the action of a generated soliton looks similar to that 
described in the previous sections for ideal solitons: in 
particular, for T = 0 an almost complete flipping is ob¬ 
tained. More pronounced differences emerge for T 7 ^ 0, 
where the asymptotic value n^(rf) is no more constant in 
time but fluctuates, being subjected to the thermal fluc¬ 
tuations of the spin chain. However, we note that such 
fluctuations are conceptually different from the decoher¬ 
ence phenomena commonly met when dealing with open 
quantum systems, as the (whatever noisy) effective field 
acting on the qubit is still classical, keeping the qubit 
evolution on the Bloch sphere. 


V. CONCLUSIONS 


Using dynamical solitons as magnetic signals running 
through spin-chains is quite a promising prospect, that 
however needs an in-depth analysis in order to become 
a more solid possibility. In fact, besides checking intu¬ 
itions, quantitative conditions must also emerge and be 
tested, with specific reference to the realization one has in 
mind. To this respect, the scheme presented in this work 
might find several different applications, as spin chains 
are versatile models that can be used for describing the 
most diverse real situations. Let us therefore end this 
work by briefly commenting upon the conditions identi¬ 
fiable as essential in our analysis, in the case of an im¬ 
plementation based on solid-state system^^lHlll 

First of all we have numerically demonstrated that 
systems of interacting magnetic moments on one¬ 
dimensional lattices, possibly of hnite length, support 
dynamical configurations which are the discrete counter¬ 
parts of /3-solitons if A /3 » d, i.e., ^/jH/JS sin/3 1, to 
be confronted with iL 7 ^ 0 for getting u > 0. Given the 
values of S and J tipically observed in magnetic com¬ 
pounds, it is JS^ ~ 1 -7 10^ K meaning that, as = 
0.67 K/Tesla, only very large helds could break the above 
inequality, and the continuum approximation is there¬ 
fore most often justihed. Further notice that our scheme 
might also be considered in the case of dynamical soli¬ 
tons that do not require H ^ 0 for being supported by the 
Heisenberg chairP^. Moreover, we know that solitons ex¬ 
ist and run also in anisotropic spin chainPPl, which makes 
our scheme potentially efficient in the case of anisotropic 
quasi one-dimensional real compounds. In fact, this is 
quite a relevant feature when thinking of implementa¬ 
tions based on one-dimensional monatomic metal chains 
deposited on surfaceP^ where the system’s geometry in- 
evitably makes the intra-chain exchange anisotropicP^I. 

As for the issue of how dynamical solitons can be gener¬ 
ated, we have preliminary indications that our proposal 
keeps being effective when the action of the field pulse 
is not punctual, as understood by Eq. ( [I^ , but 
rather extends to a finite part of the chain end (shorter 
than the induced soliton), as required in a realistic setup; 
hence, the effect of is amplified by the effective 

number of spins it affects and its intensity can well be 
comparable with H. 

Having shown that a controlled action on the qubit 
can actually be obtained by its interaction with the 
nearby running magnetic soliton, we notice that the con¬ 
dition required by the continuum appr oximatio n is fully 
consistent with the small values of y^^H/JS that are 
found to produce a permanent variation of the qubit 
state, according to the analysis presented in Sec. |IV| 
Moreover, the energy exchanged between qubit and 
chain in the case of complete flipping (obtained by, say, 
'jH/JS = 0.05, a = 0, tan/3 = 0.2, 11 = 6=1) amounts to 
SE = h{gS+jaH) ~ 10“^ as the soliton energy is of 
the order of JS^, the chain dynamics results unaffected 
by the evolution of the qubit, essentially validating the 
‘no back-action’ approximation mentioned in Sec. |IV[ As 
for the limits dictated by the typical coherence times at¬ 
tainable in solid state qubits realizations, an additional 
relevant quantity is the time tprop required by the soli¬ 
ton to reach the qubit after its injection: for a time scale 
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FIG. 6. Contour plots of the asymptotic magnetization n^(rf) as a function of the parameters tan/3 and fj., for a = 3 and 5 = 1. 
In each panel a different value of h is considered. 


of {JS)~^ ^ 10“^^s, we can estimate tprop ~ Ins, if 
the qubit lies around 10^ lattice constants away from the 
chain end. 

Finally, it is worth noticing that the magnetic solitons 
propagation we have studied is an energy-conservative 
phenomenon and has proved to be robust against ther¬ 
mal noise up to a reduced temperature T 0.01: this 
suggests that, besides the specific proposal presented in 
this work, using solitons for transferring either classi¬ 
cal or quantum information in solid-state devices might 
strongly alleviate the heat dissipation requirements that 
seriously affect more conventional solutions, without re¬ 


quiring a highly demanding lowering of the operating 
temperature. 
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